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We study strain localization as an enhanced velocity 
weakening mechanism on earthquake faults. Fault friction is 
modeled using Shear Transformation Zone (STZ) Theory, a 
microscopic physical model for non-affine rearrangements in 
granular fault gouge. STZ Theory is implemented in spring 
slider and dynamic rupture models of faults. We compare 
dynamic shear localization to deformation that is uniform 
throughout the gouge layer, and find that localized slip en- 
hances the velocity weakening of the gouge. Localized elas- 
todynamic ruptures have larger stress drops and higher peak 
slip rates than ruptures with homogeneous strain. 



1. Introduction 

The earthquake problem spans a vast range of length 
and time scales, from contacts between individual grains 
up through tectonic networks of faults. One of the primary 
modeling challenges involves identifying the relevant physi- 
cal instabilities, and accurately and efficiently propagating 
this information between scales. In this study, we investi- 
gate gouge-scale strain localization as a physical mechanism 
for enhanced velocity weakening and apply it to the larger 
scale problem of dynamic rupture propagation on faults. 
Strain localization has been observed in a variety of con- 
texts, including numerical simulations of gouge [Morgan and 
Boettcher, 1999], experimental studies of laboratory faults 
[Marone, 1998], and field observations of faults [Chester and 
Chester, 1998]. 

The formation of shear bands is a dynamic instability that 
we resolve numerically by incorporating a gouge layer of fi- 
nite width. This approach differs from the common practice 
of implementing a slip weakening or rate and state friction 
law on a planar fault. In our model, the dynamic behav- 
ior of the friction law determines how strain is distributed 
throughout the gouge layer. Within the gouge, the material 
is governed by Shear Transformation Zone (STZ) Theory, 
a microscopic physical model for non-affine deformation in 
amorphous materials [Falk and hanger, 1998, 2000]. We 
study the formation of localized shear bands, in contrast 
with homogeneous deformation. For homogeneous deforma- 
tion, strain is uniform throughout the gouge layer. In shear 
bands, slip spontaneously localizes to an interface which is 
narrow even on the scale of the gouge. The dynamic instabil- 
ity associated with localization enhances the velocity weak- 
ening that occurs with homogeneous deformation. We find 
that the enhanced weakening due to shear strain localiza- 
tion decreases the dynamic sliding friction and increases the 
peak slip rates in dynamic earthquake models. The results 
illustrate that gouge-scale strain localization has fault-scale 
consequences. 

2. STZ Friction Law 

STZ Theory is a continuum model for amorphous ma- 
terials, suitable for predicting the larger scale constitutive 
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behavior of a layer of fault gouge. STZ Theory captures 
certain features of a wide array of amorphous materials, in- 
cluding deformation of metallic glasses [Falk and hanger, 
1998, 2000], boundary lubrication [hemaitre and Carlson, 
2004], granular flow [hois et ai, 2005], shear bands in glassy 
systems [Manning et ai, 2008], and frictional weakening in 
dynamic earthquake rupture [Daub and Carlson, 2008]. 

A schematic of the fault is shown at the left in Fig. 1. 
A layer of fault gouge, which is modeled using STZ Theory, 
is sheared between elastic rocks. The plastic strain rate 7 P ; 
can vary in the z-direction in the gouge (center picture in 
Fig. 1), though the shear stress r is assumed to be uniform 
in the layer. 

Non-affine rearrangements in the gouge occur in localized 
regions, called shear transformation zones (STZs). STZs are 
modeled as bistable zones that switch between two orienta- 
tions under shear stress. An STZ can only flip once in the 
direction of the applied shear stress, so STZs are continually 
created and annihilated to sustain plastic flow. A schematic 
of an STZ switching orientation is shown at the right in 
Fig. 1. 

The plastic strain rate is influenced by two properties of 
the gouge: the number of STZs, and how these STZs change 
orientation. The number of STZs follows a Boltzmann dis- 
tribution, with an effective disorder temperature x [hanger, 
2004], and the STZs also change orientation due to the ap- 
plied stress. These two contributions can be summarized 
as: 



7 P i =exp(-l/x) R(t). 



(1) 



The function R(r) describes how the rate at which STZs 
change orientation depends on the applied shear stress. The 
physics behind our choice for R(r) is an Eyring model 
[Eyring, 1936], which matches the logarithimic velocity de- 
pendence of laboratory faults. However, the enhancement 
of velocity weakening that strain localization produces is in- 
dependent of the details of R(t). 

Regions with higher effective temperature have more 
STZs, and undergo more plastic strain. The effective tem- 
perature is distinct from the thermal temperature, but is 
similar to the free volume [hanger, 2007], which was used as 
the state variable in previous STZ friction models [hemaitre, 
2002]. However, we expect the effective temperature to ex- 
hibit dynamic behavior similar to thermal temperature, and 
we therefore include shear heating and diffusion terms in its 
governing partial differential equation. The friction equa- 
tions and the parameters for our simulations can all be found 
in Table 1. The details of the derivation can be found in 
Manning et al. [2007] and hemaitre [2002]. 

3. Spring Slider Model 

To investigate the dynamics of localization, we study the 
STZ friction law (Table 1) coupled to a non-inertial single 
degree of freedom elastic slider. A layer of gouge of width 
2w separates two blocks, one of which is pulled by a spring 
of stiffness k at a velocity Vq. We numerically solve for the 
stress and the effective temperature in the gouge layer using 
finite differences with an explicit two step time integration 
scheme. Boundary conditions on the effective temperature 
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are periodic, as the effective temperature dynamics are not 
particularly sensitive to the boundary conditions. We drive 
the slider from rest to a seismic slip rate of Vq = 1 m/s with 
a spring of stiffness k = 100000 MPa/m. This approximates 
the rapid loading and slip acceleration that occurs during 
dynamic seismic slip. We set the half width of the gouge 
layer as w = 1 m and start with an initial shear stress of 
r(t = 0) = 70 MPa. 

A special case of the spring slider dynamics arises for ho- 
mogeneous initial conditions for the effective temperature. 
In this case, by symmetry all subsequent deformation is ho- 
mogeneous. We set the initial effective temperature to be 
uniformly x(i = 0) = 0.018. A plot of stress as a function 
of displacement is illustrated in Fig. 2, which shows that 
stress weakens with displacement in a manner similar to 
the laboratory-based Dieterich-Ruina friction law [Dieterich, 
1979], and the STZ Free Volume law considered by Daub and 
Carlson [2008]. 

In general, we expect the initial conditions for effective 
temperature in the gouge to reflect the inherent heterogene- 
ity of fault zones [Chester and Chester, 1998]. A localized 
shear band dynamically forms in this case, due to the dy- 
namic response of the friction equations. A shear band, with 
a width that scales with the diffusion constant V^D, develops 
and ultimately accomodates all of the strain in the material. 

While a narrow shear band results from any non-uniform 
initial conditions, we focus on an idealized scenario to il- 
lustrate the dynamic evolution of the shear band. We add 
a small amplitude (Ax = 2 x 10~ 7 ) symmetric step per- 
turbation in the center of the gouge layer of half width 
0.1 m to the otherwise uniform initial effective temperature 
( X (t = 0) = 0.018). 

We compare the shear stress as a function of displace- 
ment for both dynamically localized strain and homoge- 
neous strain in Fig. 2(a). Uniform strain requires about 
0.3 m of displacement for the shear stress to stabilize to 
t — 63.38 MPa. Dynamic localization drops the sliding 
stress to r = 62.14 MPa much more rapidly, and illus- 
trates that localization enhances the velocity weakening of 
the gouge. 

Figures 2(b) and 2(c) show snapshots of the strain rate in 
the gouge layer for dynamically localizing strain at a series 
of representative points along the stress versus displacement 
curve. During the earliest stages of displacement, strain oc- 
curs nearly uniformly throughout the gouge layer (plot (1) 
in Fig. 2(b)). The effect of the initial perturbation is neg- 
ligible, and the strain rate is uniform. The displacement 
before localized strain begins depends on the magnitude of 
the perturbation in the initial effective temperature. This 



Table 1. Equations and parameters for the STZ constitutive 
law. The equations include the specific version of Eq. (1) and 
the partial differential equation governing the effective tem- 
perature. 



Equations 


■y p l = exp (-1/x) 2e/t exp (-/ ) (1 - t v /t) cosh (r/a d ) 
dx/dt = %iT/(c Q T y ) (l - xlog(WV)/A>) + %i D d 2 x/dz 2 


Parameter 


Description 


e = 10 
t = 10~ 6 s 

fo = 118 
Ty = 50 MPa 
u d = 0.5 MPa 

CO = 1 
70 = 80000 s" 1 

Xw = 0.5 
D = 10~ 5 m 2 


Typical number of particles in an STZ 

STZ rearrangement time scale 
STZ rearrangement activation energy 
Yield stress (below t v , 7 p ; = 0) 
Stress fluctuation for STZ rearrangement 
Effective temperature specific heat 
Strain rate at which x diverges 
Stress weakens with strain rate if Xw < 1 
Diffusion constant 



is because a region with an elevated effective temperature 
also has a higher strain rate (Eq. (1)). The shear heating 
term in the effective temperature equation (Table 1) is pro- 
portional to the strain rate, so regions of increased effective 
temperature heat more rapidly. Therefore, for larger ini- 
tial perturbations, the local effective temperature increases 
more rapidly, and less displacement is needed before the 
strain dynamically localizes. 

The strain rate profile once localization begins is shown 
in curve (2) in Fig. 2(b). The dynamic instability in the 
friction law causes strain to localize to a shear band, and 
friction weakens rapidly. A narrower shear band (plots (3) 
in Fig. 2(b) and (4) in Fig. 2(c)) develops with further dis- 
placement. By symmetry, the narrow shear band occurs in 
the center of the gouge, and its width is determined by the 
diffusion constant. Once the strain has completely localized 
to this diffusion-limited shear band, the frictional strength 
stabilizes. While the shear heating and diffusion terms do 
not exactly cancel each other, there is no noticeable change 
in the effective temperature on the time scale of an earth- 
quake rupture and the gouge slides stably. 

4. Dynamic Ruptures 

In this section, we investigate the impact of strain lo- 
calization on the propagation of spontaneous elastodynamic 
ruptures. We model the fault gouge as a thin layer with half- 
width w between two homogeneous, isotropic, linear elastic 
solids. We solve for the elastodynamic response for 2D anti- 
plane slip using a spectral boundary integral method [Perrin 
et at, 1995]. The elastodynamic equation is solved simulta- 
neously with the friction law (Table 1). 

The spectral boundary integral method can only accomo- 
date explicit time steps. The diffusion term in the effective 
temperature equation imposes a time step restriction that 
makes implementation in a rupture code impractical. How- 
ever, ignoring diffusion in the spring slider model leads to 
nearly the same stress response as a function of displace- 
ment. Therefore, as a first effort to incorporate localization 
into dynamic rupture models, we omit the diffusion term in 
the effective temperature evolution equation. 

We consider a simple fault 2 km in length along strike, 
with a gouge half width of 1 m. The friction parameters in 
Table 1 are spatially uniform along strike and throughout 
the gouge width (not including the effective temperature, 
which we solve for dynamically). The shear stress on the 
fault is initially r(i = 0) = 70 MPa, except for a small 
patch of width 0.2 km at the center of the fault where the 
stress is 79 MPa to initiate rupture. We consider two dif- 
ferent rupture scenarios, analogous to those in the spring 
slider section: one where shear strain localizes dynamically, 
and one where deformation is homogeneous. The initial ef- 
fective temperature does not vary along strike, and within 
the gouge layer it is the same as for the corresponding spring 
slider model for each rupture scenario. 

A plot comparing how stress weakens with slip for the 
two ruptures is shown in Fig. 3(a). During the initial 
stages of slip, the curves are indistinguishable. For the later 
stages of slip, shear stress weakens much more rapidly due to 
the dynamic instability of localization, increasing the stress 
drop. Homogeneous deformation and dynamic localization 
produce very different slip rates, as can be seen in Fig. 3(b). 
The rupture front arrives earlier and has a higher peak slip 
rate when strain localization occurs. 

To illustrate the importance of the dynamic instability, 
we contrast our results for localized shear with two addi- 
tional homogeneous ruptures of different fixed gouge widths 
w chosen to match particular aspects of the localized rup- 
ture. Slip velocity as a function of time is plotted at a point 
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0.35 km from the hypo-center for all models in Fig. 3(b). 
The properties that we compare are the peak slip veloc- 
ity and the time at which slip initiates. The rupture front 
in the intermediate model (w = 0.375 m) matches the ar- 
rival time of the localized rupture, but the peak slip rate 
is smaller. For the narrowest gouge thickness (w = 0.1 m), 
we see peak slip rates similar to the localized rupture but 
earlier arrival. Figure 3(c) plots the slip rate as a function 
of slip for the homogeneous rupture with w = 0.1 m and the 
localized rupture. This clearly shows that the initial broad 
deformation in the localized rupture does not simply delay 
the rupture, but also lessens the slip acceleration during the 
earliest stages of slip. 

5. Discussion 

The main effect of dynamic strain localization in our sim- 
ulations is to provide a mechanism for enhanced velocity 
weakening. The dynamic instability that produces localiza- 
tion has two effects on the shear stress: a reduced dynamic 
sliding stress, and a very rapid decrease in stress with slip. 
Allowing for enhanced weakening makes our faults more un- 
stable - both the stress drop and the peak slip rates increase 
for dynamic localization relative to homogeneous slip. 

The width of the narrow shear band in our simulations is 
determined by the diffusion length scale \/Z). For different 
driving rates, other shear band widths are possible. 

Other weakening mechanisms have been discussed in the 
context of rapid seismic slip, including pore fluid prcssuriza- 
tion [Lachenbruch, 1980], flash heating at asperity contacts 
[Tullis and Goldsby, 2003], production of a thixotropic silica 
gel [Di Toro et al, 2004], and frictional melting [Di Toro et 
al., 2006]. Most of these are due to (regular) thermal effects. 
Localization complements these, as localization occurs be- 
fore significant change in the thermal temperature. Rupture 
models that employ both could examine how thermal weak- 
ening and shear localization interact. 

The enhanced weakening effect of dynamic localization 
occurs regardless of the choice of the function R(t) in 
Eq. (1). Other forms of R(t) will still exhibit more rapid 
weakening during localized deformation. This indicates that 
shear localization can affect a variety of materials, and have 
consequences for other geophysical systems. 
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Figure 1. An illustration of the range of scales in the 
earthquake rupture problem. The scale moves to progres- 
sively smaller phenomena from left to right. Left: fault 
scale model, with a thin layer of gouge described by STZ 
theory sheared between elastic rocks. Center: close up of 
deformation inside the gouge, where shear strain develops 
into a localized shear band (dark regions). Shear band 
image taken from Folk and Shi [2002], and re-oriented to 
match the sense of shear of the fault and grains. Right: 
microscopic picture of the grain scale, with an STZ un- 
dergoing transformation from a "+" oriented zone (left) 
to a "— " oriented zone (right). As the gouge deforms 
plastically, the ellipse drawn through the particles flips 
its orientation. STZ diagram taken from Falk and hanger 
[1998]. 
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Figure 2. (a) Plot of shear stress as a function of shear 
displacement for the spring slider model. Comparison 
between dynamic strain localization and homogeneous 
deformation reveals that localized strain exhibits more 
rapid weakening, (b) Strain rate profiles for four rep- 
resentative shear displacements during localized strain, 
indicated in plot (a). (1) An initial period of broad defor- 
mation occurs before (2) strain dynamically localizes. (3) 
A narrower, diffusion-limited shear band develops with 
further shear, (c) Inset: The narrow shear band eventu- 
ally accommodates all of the deformation. 
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Figure 3. Dynamic rupture evolution at a point 0.35 km 
from the hypocenter. (a) Comparison of shear stress as 
a function of slip. Dynamic localization of deformation 
produces more rapid velocity weakening than the rup- 
ture with homogeneous strain, (b) Plot of slip rate as 
a function of time. The dynamic strain localization rup- 
ture is compared with a host of models with homogeneous 
strain. None of the values of the imposed gouge width w 
can match both the peak slip rate and rupture front ar- 
rival of the rupture with localized strain, (c) Inset: Slip 
rate as a function of slip for the localized and narrow- 
est width homogeneous rupture. The more rapid accel- 
eration of slip in the narrowest homogeneous rupture is 
distinct from the localized model. 



